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Abstract 

The leading non-perturbative contribution to the static QCD potential at r -C Aq^j-, is 
known to be O(r^) in operator-product expansion. It indicates that a "Coulomb+linear" 
potential at r ^ ^qcd included in the perturbative QCD prediction of the potential. 
It was shown that this is indeed the case, and the "Coulomb+linear" potential has been 
systematically computed up to next-to-next-to-leading logarithmic order, in terms of 
the Lambda parameter in the MS scheme (Aj^). We review the present status of the 
perturbative prediction for the QCD potential, which takes into account the theoretical 
breakthrough that took place around 1998. 



*Talk given at the "13th International Seminar on High Energy Physics (Quarks 2004)," Pushkinskie 
Gory, Russia, May 24 - 30, 2004. 



1 Introduction 



In this article we review the present status of perturbative QCD predictions for the QCD 
potential, in the distance region relevant to the bottomonium and charmonium states, 
0.5 GeV~^(0.1 fm) ^ r < 5 GeV~^(l fm). We take into account the theoretical breakthrough 
that took place around 1998 [1], which improved the accuracy of perturbative prediction dra- 
matically. 

For a quite long time, the perturbative QCD predictions of the QCD potential VqcdI'^) 
were not successful in the above distance region. In fact, the perturbative series turned out 
to be very poorly convergent at r > 0.5 GeV~^; uncertainty of the series is so large that one 
could hardly obtain meaningful predictions. Even if one tries to improve the perturbation 
series by certain resummation prescription (such as renormalization group improvement), 
scheme dependence of the result was also very large [2]; hence, one could neither obtain 
accurate prediction of the potential in this distance region. It was later pointed out that the 
large uncertainty of the perturbative QCD prediction can be understood as caused by the 
leading-order [C(Aqcd)] infrared (IR) renormalon contained in VqcdI'^) [3]. 

On the other hand, empirically it has been known that phenomenological potentials 
and lattice computations of Vqcd('") are both approximated well by the sum of a Coulomb 
potential and a linear potential in the above range 0.5 GeV-^ < r < 5 GeV'^ (See e.g. 
[4]). The linear behavior of Vqcd('") at large distances r ^ ^qcd (verified numerically by 
lattice computations) is consistent with the quark confinement picture. For this reason, and 
given the very poor predictability of perturbative QCD, it was often said that, while the 
"Coulomb" part of Vqcd('") (with logarithmic correction at short-distances) is contained in 
the perturbative QCD prediction, the linear part is purely non-perturbative and absent in 
the perturbative QCD prediction even at r < Aqcd' t\\.aX the linear potential needs 
to be added to the perturbative prediciton to obtain the full QCD potential. See Fig. 1. 
Nevertheless, to the best of our knowledge, there was no firm theoretical basis for this 
argument. 

Several years ago, the perturbative prediction of the QCD potential became much more 
accurate in this region. There were two important developments: (1) The complete 0{a\) 
corrections to the QCD potential have been computed [5]; also the relation between the 
quark pole mass and the MS mass has been computed up to 0(tt|) [6]. (2) A renormalon 
cancellation was discovered [1] in the total energy of a static QQ pair*, -E'tot(^) = 2mpoie + 
Vqcd('^)- Consequently, convergence of the perturbative series of Ex,ot{j) improves drastically, 
if it is expressed in terms of the quark MS mass instead of the pole mass. 

Physically, improvement of convergence stems from decoupling of IR gluons from the 
color-singlet QQ system. Intuitively, IR gluons, whose wavelengths are larger than r, cannot 
see the color charge of Q or Q but only see the total charge of this system. Therefore, in the 

*By "static", we mean that the kinetic energies of Q and Q are neglected. 
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Figure 1: Status of theoretical predictions for the QCD potential before the discovery of 
renormalon cancellation. One could hardly see sign of convergence of perturbative series at 



computation of the total energy of the QQ system, we expect that contributions of IR gluons 
should decouple. This is naturally realized, if we use a quark mass, which is renormalized 
to include contributions of only UV gluons, such as the MS mass; in this case, there is a 
cancellation of IR contributions between the self-energies oi Q/Q and the potential energy. 
If we use a mass, which includes also contributions of IR gluons, such as the pole mass,^ 
then the decoupling of IR gluons is not realized. Generally, convergence of a perturbative 
expansion is worse when there are more contributions from IR gluons, and vice versa. 

Let us demonstrate the improvement of accuracy of the perturbative prediction for the 
total energy Etot{r) up to 0{ag), when the cancellation of C(Aqcd) renormalons is incor- 
porated. As an example, we take the bottomonium case: We choose the MS mass of the 
6-quark, renormalized at the 6-quark MS mass, as frih = m^^(m^^) = 4.190 GeV; in internal 
loops, four flavors of light quarks are included with m„ = fnd = rris = and nic = 1.243 GeV. 
(See the formula for E^^^{r) in [7].) In Fig. 2, we fix r = 2.5 GeV~^ (midst in the distance 
range of our interest) and examine the renormalization-scale (/z) dependences of Etot{f)- We 
see that Etot{r) is much less scale dependent when we use the MS mass than when we use the 
pole mass. Hence, the perturbative prediction of Etot{i^) is much more stable in the former 
scheme. 

tPole mass is defined as the energy of a free quark at rest. Hence. IR gluons can sec its color charge and 
contribute to the pole mass. 
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Figure 2: Scale dependences of E^^^{r) up to 0(a|) at r = 2.5 GeV~^, in the pole-mass and 
MS-mass schemes. A horizontal line at 8 GeV is plotted as a guide. 



We also compare the convergence behaviors of the perturbative series of E^ot for the 
same r and when /i is fixed to the value, at which E'tot becomes least sensitive to variation 
of /i (minimal-sensitivity scale). Convergence of the perturbation series turns out to be 
close to optimal for this scale choice: at r = 2.5 GeV~^, the minimal-sensitivity scale is 
/i = 0.90 GeV. We find 

E^l^{r) = 10.408 - 0.275 - 0.362 - 0.784 GeV (Pole-mass scheme) (1) 
= 8.380 + 1.560 - 0.116 - 0.022 GeV (MS-mass scheme). (2) 

The four numbers represent the 0{ag), 0{ag), 0(a|) and C?(a|) terms of the series expan- 
sion in each scheme. The 0{a^g) terms represent merely the twice of the pole mass and of 
the MS mass, respectively. As can be seen, if we use the pole mass, the series is not con- 
verging beyond 0{ag). On the other hand, in the MS-mass scheme, the series is converging. 
One may also verify that, when the series is converging (MS-mass scheme), /i-dependence of 
-E'tot('") decreases as we include more terms of the perturbative series, whereas when the se- 
ries is diverging (pole- mass scheme), /^-dependence does not decrease with increasing order. 
(See e.g. [8].) 

We observe qualitatively the same features at different r and for different number of light 
quark flavors n/, or when we change values of the masses ffib, fUc- Generally, at smaller r, 
Etot{r) becomes less /i-dependent and more convergent, due to asymptotic freedom. See [7,9] 
for details. 

The aim of this paper is to review the properties of Etot{r), given the much more ac- 
curate prediction as compared to several years ago. In Sec. 2 we examine Etot{f) up to 



3 




r [GcV^M 



Figure 3: Etot{r) and typical phenomenological potentials. Constants are added to make all curves 
coincide at r = 1 GeV~^. Arrows at the bottom show the r.m.s. radii of the heavy quarkonium 
states. The figure is taken from [7]. 

0{ag). Sec. 3 provides a classification of perturbative and non-perturbative contributions 
to Etot{r) in terms of operator-product-expansion (OPE). We present our main result, that 
the perturbative prediction of the leading short-distance contribution to Etotif) is given as 
a "Coulomb+linear" potential, in Sec. 4. Conclusions are given in Sec. 5. 

2 Pert, prediction of Etot{r) up to 0(a|) 

Let us first review comparisons of the perturbative predictions of Etot{r) up to with 
phenomenological potentials and with lattice computations of the QCD potential. 

In Fig. 3, £'tot('") is compared with typical phenomenological potentials. Since the latter 
are determined from the heavy quarkonium spectra, realistic values for the input parameters 
of Etoti'r), as given in the previous section, have been chosen. The scale /i = /i(r) is fixed by 
either dEtot/dfi = (minimum sensitivity scale) or \El^l\ = minimum, but both prescrip- 
tions lead to almost same values of Etot{i^)- See [10,7,9] for details. We see that Etot{r) 
corresponding to the present values of the strong coupling constant (dashed lines) agree well 
with the phenomenological potentials within estimated perturbative uncertainties (indicated 
by error bars). We also note that the agreement is lost quickly if we take as{Mz) outside of 
the present world-average values, so that the agreement is unlikely to be accidental. 

In fact, by now several works have confirmed agreements of the perturbative predictions of 
the QCD potential with phenomenological potentials or with lattice computations [10,11,7, 
12,13,9]. Although details depend on how the renormalon in the QCD potential is cancelled. 
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qualitatively the same conclusions were drawn: i.e. perturbative predictions become accurate 
and agree with phenomenological potentials/lattice results up to much larger r than before. 
In particular, in the differences between the perturbative predictions and phenomenological 
potentials/lattice results, a linear potential of order AqQ^r at distances r < Aq^j-, was ruled 
out numerically. In other words, one cannot get the full QCD potential if one adds a linear 
potential to the perturbative QCD potential (after renormalon cancellation). 

A crucial point is that, once the (9(Aqcd) renormalon is cancelled and the perturbative 
prediction is made accurate, the perturbative potential becomes steeper than the Coulomb 
potential as r increases. This feature is understood, within perturbative QCD, as an effect 
of the running of the strong coupling constant [14,10,11]. In short, the interquark force 
becomes stronger than the Coulomb force at larger r by this effect. 



3 Operator-product expansion of Vqcd(^) 

Operator-product expansion (OPE) of Vqcd('") for r <^ Aq^^ was developed [15] within an 
effective field theory "potential non-relativistic QCD" (pNRQCD) [16]. In this framework, 
Vqcd('") is expanded in r (multipole expansion), when the following hierarchy of scales exists: 

1 

AQCD</i/<-. (3) 

Here, /ij denotes the factorization scale. At each order of the expansion in r, short-distance 
contributions {q > Hf) are factorized into perturbatively computable Wilson coefficients and 
long-distance contributions (g < ///) into matrix elements of operators, i.e. non-perturbative 
quantities. The leading non-perturbative contribution to the potential turns out to be 

Explicitly, the QCD potential is given by 



VqcDir) = Vs{r; fij) + 6Eus{r; ///), (4) 

rji poo 

^^c Jo 

+ Oir% (5) 

The leading short-distance contribution to Vqcd('") is given by the singlet potential Vs{r). 
It is a Wilson coefficient, which represents the potential between the static QQ pair in color 
singlet state. The leading non-perturbative contribution [O(r^) in the multipole expansion] 
is contained in the matrix element in eq. (5). AV{r) = Vo{r) — Vs{r) denotes the difference 
between the octet and singlet potentials; see [15] for details. 

The singlet potential Vs{r]fif) can be computed in perturbative expansion in as by 
matching pNRQCD to QCD. It is essentially the perturbative expansion of Vqcd('") after 
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the contribution of soft degrees of freedom (g < fif) is subtracted. In particular, the IR 
divergences and IR renormalons contained in the perturbative expansion of VqcDif) have 
been subtracted via appropriate renormalization prescription [15,12,18]. Thus, Vs{r] fif) can 
be computed accurately by perturbative QCD. 

We may understand why the leading non-perturbative contribution is (9(r^) as follows. 
As well known, the leading interaction (in expansion of r) between soft gluons and a color- 
singlet QQ state of size r is given by the dipole interaction r ■ E"-, where E°- denotes the color 
electric field. It turns the color singlet QQ state into a color octet QQ state by emission of 
soft gluon(s). To return to the color singlet QQ state, it needs to reabsorb the soft gluon(s), 
which requires an additional dipole interaction. Thus, the leading contribution of soft gluons 
to the total energy is 0{r^). 



4 Pert, prediction of Vs{r) up to NNLL: 
"Coulomb+linear" potential 

Here, we present our main result. As already stated, the singlet potential Vs{r;fj,f) can be 
computed in perturbative QCD, which may be improved via renormalization group (RG). We 
have shown that Vs{r; fif) thus computed can be expressed as a "Coulomb+linear" potential, 
up to an (9(r^) correction. The correction can be absorbed into the (9(r^) non-perturbative 
contribution SEus{r; fif) by appropriate renormalization prescription. Explicitly, we find 
[17,18] 

Vs{r; fif) = const. + Vdr) + ar + 0{fi}r^), (6) 

where 

Vc{r) = Im / dq (7) 

a = ^ [ dqqaviq). (8) 

'C2 



2111 



In the above equations, av{q) denotes the perturbative evaluation of the l^-scheme coupling 
constant in momentum space [5] (improved by RG, after subtraction of IR divergences via 
renormalization). The integral paths Ci and C2 are displayed in Figs. 4. Vc(r) and a are 
defined to be independent of fif (at least) up to NNLL. Cp = 4/3 denotes a color factor; 
/3o = 11 — 2?7,i/3 represents the 1-loop coefficient of the beta function. In eq. (6), we are not 
concerned about the constant (r-independent) part of Vs{r), keeping in mind that it can 
always be absorbed into 2mpoic in the total energy E\_q^_{t). 
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Figure 4: Integral contours C\ and C2 on the complex g'-plane. denotes the Landau singularity 
of For 1-loop running, is a pole; for 2- and 3- loop running, g,, is a branch point. In the 

latter case, branch cut is on the real axis starting from q^^ to —00. 



By evaluating the above integral, the coefficient of the linear potential a can be expressed 
analytically in terms of the Lambda parameter in the MS-scheme A^. For instance, at LL, 
it is given by 

2tiCf ( . l_loop\^ 



The "Coulomb" potential has a short-distance asymptotic behavior consistent with RG, 

yc(r) ~ ^"f" . (r ^ 0). (10) 



whereas its long-distance behavior is given by 

Vc{r) ~ (r ^ 00). (11) 

In the intermediate region both asymptotic forms are smoothly interpolated. 

Fig. 5 shows the Coulomb+linear potentials corresponding to the RG improvement by 
the 1-loop running (LL), 2-loop running with 1-loop non- logarithmic term (NLL), and 3- 
loop running with 2-loop non- logarithmic term (NNLL). They are compared with the lattice 
results. Since the 3-loop non-logarithmic term is not yet known, the 3-loop running case 
represents our current best knowledge. Up to this order, the Coulomb-|-linear potential 
agrees with the lattice results up to larger r as we increase the order.* 

*The NNLL originating from the ultra-soft scale [19] hardly changes the 3-loop running case displayed in 
Fig. 5 [17]. 
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Figure 5: Comparison of the Coulomb+lincar potential (n; ~ 0) with the lattice results in the quenched 
approximation. 1-, 2- and 3-loop correspond, respectively, to LL, NLL and NNLL in the RG-improved 
computation of Vs (r) . 

5 Conclusions 

After discovery of renormalon cancellation, perturbative prediction of the total energy £'tot('") = 
2?Tipoie+VQCD('^) became much more accurate in the distance range 0.5 GeV~^ ^ t < 5 GeV~^. 
Consequently we observe the following: 

• Perturbative prediction of Etot{r) up to agrees well with phenomenological 
potentials/lattice results within the estimated perturbative uncertainty. 

• According to OPE, the leading non-perturbative contribution to Etot{r) in expansion 
in r is (9(r^). 

• Perturbative QCD prediciton of the singlet potential Vs{r) (leading Wilson coeffi- 
cient in OPE) can be written as a "Coulomb+linear" potential. Up to NNLL, the 
Coulomb+linear potential shows a converge towards lattice results. 

All the theoretical analyses conclude that, there is no freedom to add a linear potential of 
0(AQQj-,r) to the perturbative prediction of the QCD potential at r < Aq^q. Therefore, if 
we are able to define the string tension from the shape of the potential at r ^ AqJ^j-, (it has 
been empirically the case in phenomenological potential model approach), the string tension 
may be within the reach of perturbative prediction. 
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